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FOREWORD 


The  objective  of  the  program  reported  herein  was  to 
provide  the  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base,  Ohio,  with  a  background  on  Monte  Carlo  time 
domain  simulation  techniques  for  application  to  aerospace 
vehicle  structures.  Dr.  Jon  Lee  of  WRDC/FIBG  was  the  con¬ 
tract  monitor. 

The  work  reported  was  performed  for  the  Aerospace 
Structures  Information  and  Analysis  Center  (ASIAC),  which 
is  operated  for  the  Flight  Dynamics  Laboratory,  Wright 
Research  and  Development  Center,  Wright-Patterson  Air  Force 
Base,  Ohio,  by  Anamet  Laboratories,  Inc.  under  Contract  No. 
F33615-87-C-3228.  The  ASIAC  Task  Number  was  4.2-34. 
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1.0  INTRODUCTION 


The  prediction  of  dynamic  response  and  fatigue  life  of  surface  protec¬ 
tion  systems  of  high  speed  aircraft  is  a  critical  task  in  the  development  of  modem  flight 
structures.  With  an  increasing  use  of  lightweight  materials,  such  as  fiber-reir\forced 
composites,  metal-matrix  composites,  and  various  high  temperature  resistant 
intermetallic  compounds  on  new  generation  aircraft,  orbital  and  sub-orbital  vehicles, 
there  is  a  need  for  improved  methods  for  predicting  dynamic  response  and  sonic 
fatigue  life.  The  future  high  speed  supersonic/hypersonic  vehicles  will  be  subjected  to 
severe  aerodynamic  surface  environment,  intense  acoustic  loads  and  temperatures 
exceeding  3000*^  F.  The  structural  and  various  control  surfaces  are  expected  to  behave 
in  a  highly  nonlinear  fashion  with  respect  to  both  geometry  and  materials. 

A  time-domain  Monte  Carlo  approach  could  provide  the  r^uired  methodology 
for  the  solution  of  the  various  nonlinear  problems  in  relation  to  thermal  surface 
protection  systems  of  high  speed  aircraft.  These  types  of  procedures  have  been  utilized 
for  response  and  fatigue  analysis  on  a  variety  of  structural  components  exposed  to 
severe  aerodynamic,  acoustic  and  thermal  environment  [1-10].  To  implement  the  time 
domain  approach  for  nonlinear  problems  with  random  inputs,  simulation  of  random 
surface  pressures  in  time-space  domain  are  required.  A  computer  code  called 
SIMLOAD,  which  simulates  in  time  domain  a  single  dimensional  and  single  variate 
random  pressure,  has  been  implemented  on  a  WRDC  computer.  These  simulated 
time  histories  are  the  required  input  to  the  nonlinear  time  domain  problem.  One  of 
the  basic  structural  components  of  a  thermal  surface  protection  system  is  a  rectangular 
plate.  A  nonlinear  deflection/stress  response  solution  has  been  developed  in  time 
domain  for  a  rectangular  panel.  The  computer  code  called  PLATE,  which  solves  for 
nonlinear  deflection  and  nonlinear  stress  under  uniformly  distributed  random 
pressures,  was  also  implemented  on  a  WREXl  computer. 


2.0  TECHNICAL  DISCUSSION 

A  detailed  analysis  on  the  simulation  of  multidimensional/multivariate 
stationary/ Gaussian  random  processes,  solution  of  nonlinear  time  domain  random 
problems,  and  sonic  fatigue  using  time  domain  approach  has  been  presented  to 
researchers  at  WRIX2,  WPAFB.  In  what  follows,  a  brief  summary  of  these  procedures 
is  given. 


2.1  Simulation  of  Random  Input  Pressures  in  Space-Time  Domain 

Consider  a  stationary,  homogeneous  and  Gaussian  random  pressure  p(x,y,t) 
acting  on  the  surface  of  a  high  sp)eed  flight  vehicle.  The  pressure  acting  normal  to  the 
surface  varies  randomly  in  time  and  space  along  the  surface  coordinates  x  and  y.  The 
random  pressure  is  characterized  by  a  cross-spectral  density  function  SpC^.T],©),  where 
^  =  X|  -  X2  and  t]  =  yi  -  y2  are  the  spatial  separations  and  O)  is  frequency.  The  spectral 

density  of  pressiire  p(x,y,t)  can  be  obtained  in  wave-number-frequency  domain  by 
taking  the  Fourier  transformaton  of  Sp(q,ii,co)  as 


5  ^ ) ~ 
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(1) 


in  which  k|  and  k2  are  the  wave  numbers.  The  random  sample  functions  of  pressure 
p(x,y,t)  with  a  zero  mean  can  be  generated  utilizing  simulation  procedures  of 
stationary-Gaussian  random  processes  and  Fast  Fourier  Transform  (FFT)  techniques 
[11,12].  Then,  p(x,y,t)  can  be  simulated  by  the  series 
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where 


■^</r  “[2'Sp(A:ij,/fc2y,co^)A^iAA:2Aou] 


(3) 


are  realized  values  of  independent  random  phase  angles  uniformly  distributed 
between  0  and  2%,  and  Re  indicates  the  real  part.  The  values  of  the  spectra  are  selected  at 
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where  the  wave  number  and  frequency  intervals  are 


C^2u~  ^ 2 

Aa)»(a)„-cju,)/A^3 


in  which  the  subscripts  u  and  1  indicate  the  lower  and  the  upper  cut-off  values  of  the 
wave  number  and  frequency,  respectively.  The  random  pressure  p(x,y,t)  is  evaluated  at 
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=  mAx , 


/n  =  0 , 1 , . . . ,  M 1  -  1 
^  =  0, 1  ,...,^2"  1 
q*0, 1 . M3-  1 
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In  equations  (2)  and  (6) 


M,  =  2'”‘  =  v,A^i>A/i  =  2'‘' 
M2  =  2'”*  =  V2A/2>  A^2*2'** 

M3=2'”^  =  V3A/3>  A/3=2'‘’ 

Ax  =  2ji/MiA^i  =  2n/v,A:i„ 

A:/  =  2ji/M2A^2  =  2n/v2A:2„ 
Ar  =  2n/M3A(ju  =  2ji/v3a)„ 


(7) 


(8) 


If  a  random  pressure  acting  on  a  structural  surface  can  be  assumed  to  be  uniformly 
distributed  in  space,  equation  (2)  reduces  to  one-dimensional  simulation 


“M-l 
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(9) 


where 


^r  =  [25p(co^)Aou]* 


(10) 


and  Sp(co)  is  the  power  spectral  density  of  random  pressure  p(t). 

In  order  that  random  pressure  can  be  simulated  either  from  equation  (2)  or 
equation  (9),  the  input  spectral  densities  need  to  be  prescribed.  Empirical  forms  of 
cross-spectral  densities  or  spectral  densities  are  available  for  jet  noise,  turbulent 
boundary  layer  flow  pressure  and  rocket  noise,  etc.  [13-16].  Useful  approximations  can 
be  obtained  by  assuming  the  input  pressure  to  be  uniformly  distributed  over  the 
structural  component  surface.  Then,  the  cross-spectral  density  can  be  approximated  by 
band-limited  Gaussian  white  noise 


Sp{^,T],CO)  =  { 


0 


ifG)l  <  Q>  <  CO^ 

^  (o  <  CO  I  or  CO  >C0i, 


(11) 


in  which  ©l  and  ©u  indicate  the  lower  and  the  upper  cut-off  frequencies.  The 
expression  for  So  can  be  written 

(12) 

Aq) 

where  po  is  the  reference  pressure,  po  =  2.9  x  10"^  psi  (0.00002  N/m^),  A©  is  the  selected 
bandwidth,  and  SPL  is  the  sound  pressure  level  expressed  in  decibels. 
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2.2  Simulation  of  Multivariate  Random  Process 

There  are  many  problems  of  interest  where  the  time  histories  of  random  point 
forces  are  required.  Consider  a  set  of  stationary  and  Gaussian  random  processes 
Fj(t),  j  =  1,2, ...,K.  Assume  the  mean  value  is  zero  for  each  of  Fj.  Then,  the  random 
processes  can  be  simulated  as  [11,12] 


where 


Fj{t,)  =  ,/2Kd}Re\Fj{t,)\ 


(13) 


m  =  l  ;=1 


(14) 


and  the  matrix  Hjin(wi)  is  constructed  from  the  aoss-spectral  density  matrix  [Sp(a))]  of 
the  generalized  random  forces  Fj(t)  [12].  The  phase  angles  are 


(cq)  =  tarC'^ 


'  4 

Re\Hj„  (G:t)] 


(15) 


In  equation  (14),  (pml  are  the  random  phase  angles  uniformly  distributed  between  0 
and  27C, 
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Ao)  =  (OulN 
N  =  2\n  =  1,2,... 

A/  =  2'",  m  =  1,2... 

=  <jAt,  q  =  0, 1,2,..,A/-1 

Ar  =  IkKMAco)  =  2;r/|^VQji,  j 

V  =  2'"'",  m  >  n  (16) 

In  equation  (16),  N  is  the  number  of  frequency  divisions  on  cross-spectral  densities 
[Spfco)]  and  M  is  the  total  number  of  simulation  points.  The  calculation  of  matrix  {H(q))] 
can  be  found  in  Ref.  12. 


2.3  Nonlinear  Response  of  Surface  Panels  Subjected  to  High  Intensity  Noise 

The  time  domain  solutions  for  nonlinear  deflection  and  stress  response  of 
homogeneous  panels,  discretely  stiffened  panels  and  panels  made  from  composite 
materials  have  been  developed  in  Refs.  1-10.  These  solutions  were  obtained  utilizing  a 
Galerkin  type  procedure  and  numerical  integration  of  the  governing  equations  of 
motion  in  time  domain.  To  illustrate  this  procedure,  consider  a  homogeneous  panel 
shown  in  Fig.  1,  exposed  to  random  pressvire  p(x,y,t).  The  random  input  pressure  is 
assumed  to  be  of  magnitude  to  cause  the  panel  to  vibrate  in  a  nonlinear  fashion. 
Thermal,  cavity,  nonsteady  aerodynamic  and  in-plane  loads  are  not  included  in  the 
formulation.  Then,  utilizing  von  Karman’s  theory  of  large  deflection  of  thin  plates, 
the  governing  equations  of  motion  are 


dy^ 


P(x,y,t) 


(17) 
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V^F 


Eh 


d^w 

d^w' 

[[dxdy] 

dx^  ' 

dy\ 

(18) 


where  F  is  the  Airy  stress  function  and  D  =  Eh3/12  (l-v^),  E,h,  mp,v  are  modulus  of 
elasticity,  plate  thickness,  mass  per  imit  area  and  Poisson's  ratio,  respectively.  The 
membrane  in-plane  forces  Nx,Ny  and  Nxy  are  given  by 


=  d^F/dy^ 

Ny  =  d^F/dx^ 


Nxy  =  -^FldxBy 


(19) 


The  in-plane  displcements  u  and  v  are  expressed  as 


dx  Eh  ^  d^x  j 


^  _  J_  f^F  _  <92/r> 

dy  Eh  ^  B^y  j 


The  in-plane  boundary  conditions  are  satisfied  on  the  average 


a  a  du 
3  dx 


dxdy  =  0 


(20) 


(21) 


(22) 
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nb 

-dxdy  =  0 
3  dy 


1  rb 

a  ^0  ^ 


X  =0,a 


d\-  =  0 


y  =  o,b 


dx  =  0 


(23) 

(24) 

(25) 


Equations  (22)  and  (23)  imply  no  in-plane  stretching  of  the  panel  edges  (immovable 
edges)  and  equations  (24)  and  (25)  indicate  zero  shear  forces  at  the  edges. 

To  solve  equations  (17)  and  (18)  panel  normal  deflections  are  expanded  in  terms 
of  modes  6nm(x,y)  such  that 


w(x,y,t)  =  IS  ^mft  it)<l>m^(x,yy 

m  =  1  «  =  1 


(26) 


in  which  Amn  are  modal  amplitudes.  Substituting  equation  (26)  into  equation  (18)  and 
satisfying  the  in-plane  boundary  conditions,  the  solution  for  the  Airy  stress  function  F 
can  be  developed  [2-5].  Then  utilizing  this  solution,  the  Galerkin  type  procedure  and 
equations  (17)  and  (26),  a  system  of  nonlinear  differential  equations  in  the  generalized 
coordinates  Anxn  is  obtained.  After  the  simulation  of  the  random  input  pressure 
p(x,y,t),  the  nonlinear  equations  of  motion  are  solved  in  time  domain  utilizing  a 
numerical  integration  procedure. 

The  nonlinear  stresses  of  the  isothermal  panel  can  be  calculated  from 


Ez  /d^w  d^w\ 
1— 


(27) 
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Ez  /d^w  d^w\  1 

l-t/2  (28) 


E  z  d^w  Nxf 
I  +  u  dx  dy  h 


(29) 


where  z  is  the  distance  from  the  midplane  of  the  panel.  The  first  set  of  terms  in 
equations  (27-29)  correspond  to  bending  stresses  and  the  second  term  to  stresses  due  to 
midplane  stretching. 

2.4  Sonic  Fatigue 

The  nonlinear  response  leads  to  a  non-Gaussian  distribution  of  stress  amplitude. 
A  starting  point  in  developing  a  damage  theory  that  is  consistent  with  the  present  time 
domain  approach  for  nonlinear/ non-Gaussian  response  is  to  utilize  fatigue  data  from 
constant  amplitude  coupon  testing.  The  stress  versus  the  number  of  cycles  can  be 
written  as 


S  ^  =  5//^  (30) 

where  S  is  a  fixed  stress  amplitude  for  constant  amplitude  loading,  N  is  the  number  of 
stress  cycles  until  failure,  and  X  and  B  are  material  constants.  The  cumulative  damage 
can  be  written  as  [17] 


(31) 


in  which  ni  represents  the  actual  number  of  cycles  at  a  given  stress  level  and  Ni  is  the 
number  of  total  cycles  at  which  failure  occurs  at  the  same  stress  level.  Combining 
equations  (30)  and  (31) 


D=  1  £1.(5,) IS,  I" 

O  i 


(32) 


where  an  absolute  sign  is  introduced  to  account  for  positive  stress  peaks  that  might 
occur  in  the  negative  stress  region.  Let  now  D(t)  given  in  equation  (32)  represent  the 
random  accumulated  damage  per  unit  time  due  to  continuous  random  stress  S(t). 
Then,  it  can  be  shown  that  the  expected  fatigue  damage  in  a  time  interval  t  is  [18] 


E[D{x)] 


B 


\S‘{^P,{S)dS 


(33) 


where  pi(S)  is  the  probability  density  fimction  of  the  stress  peak  magnitude.  These 
probability  density  functions  are  determined  as  histograms  directly  from  the  nonlinear 
stress  response  time  history.  For  a  stationary  stress  response,  fatigue  life  can  be 
computed  from  equation  (33)  by  setting  E[D(t)]  =  1.  The  expected  total  number  of  peaks 
E[Mt]  can  be  estimated  from  [18] 


sPg^g(s  ,0 ,  s)difds 


(34) 


where  PSSS  is  the  joint  probability  density  of  S  (stress),  S(s  tress  velocity  ),  and  S  (stress 


acceleration). 
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Figure  1  A  Rectangular  Panel  Exposed  to  Random  Pressure 


